Domain wall roughness in epitaxial ferroelectric PbZro.2Tio.8O3 thin films 
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The static configuration of ferroelectric domain walls was investigated using atomic force mi- 
croscopy on epitaxial PbZro.2Tio.8O3 thin films. Measurements of domain wall roughness reveal a 
power law growth of the correlation function of relative displacements B{L) oc L^^ with ~ 0.26 at 
short length scales L, followed by an apparent saturation at large L. In the same films, the dynamic 
exponent fi was found to be ~ 0.6 from independent measurements of domain wall creep. These 
results give an effective domain wall dimensionality of d — 2.5, in good agreement with theoretical 
calculations for a two-dimensional elastic interface in the presence of random-bond disorder and 
long range dipolar interactions. 

PACS numbers: 77.80.-e, 77.80.Dj, 77.80.Fm 
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Understanding the behavior of clastic objects pinned 
by periodic or disorder potentials is of crucial importance 
for a large number of physical systems ranging from vor- 
tex lattices in type II superconductors charge den- 
sity waves and Wigner crystals to interfaces during 
growth and fluid invasion jH processes, and magnetic 
domain walls 0. Ferroelectric materials, whose switch- 
able polarization and piezoelectric and pyroelectric prop- 
erties make them particularly promising for applications 
such as non- volatile memories actuators, and sen- 

sors 10], are another such system. In these materials, 
regions with different symmetry-equivalent ground states 
characterized by a stable remanent polarization are sep- 
arated by elastic domain walls. The application of an 
electric field favors one polarization state by reducing 
the energy necessary to create a nucleus with polariza- 
tion parallel to the field, and thereby promotes domain 
wall motion. Since most of the proposed applications use 
multi-domain configurations, understanding the mecha- 
nisms that control domain wall propagation and pinning 
in fcrroelectrics is an important issue. 

A phenomenological model derived from measurements 
of domain growth in bulk fcrroelectrics (IL 12, 13] ini- 
tially suggested that the domain walls were pinned by 
the periodic potential of the crystal lattice itself. Such 
pinning was deemed possible because of the extreme thin- 
ness of ferroelectric domain walls (different from the case 
of magnetic systems}. However, measurements of the 
piezoelectric effect l22l. dielectric permittivity jl^, and 
dielectric dispersion |2J] in ferroelectric ceramics and sol- 
gel films have shown some features which cannot be de- 
scribed by the existing phenomenological theories. A mi- 
croscopic study of ferroelectric domain walls could resolve 
these issues. Recently, we have measured domain wall 
velocity in epitaxial PbZro.2Tio.8O3 thin films, showing 
that in this case commensurate lattice pining is in fact 
not the dominant mechanism 0,^1]. Rather, a creep- 
like velocity (v) response to an externally applied electric 
field E was observed with v ^ expj— C/i?^], where C is 
a constant. The exponent fi characterizing the dynamic 



behavior of the system is a function of the domain wall 
dimensionality and the nature of the pinning potential. 
These results suggested that domain wall creep in ferro- 
electric films is a disorder-controlled process. However, 
questions about the microscopic nature of the disorder 
were left open by the dynamical measurements alone. In 
order to ascertain the precise physics of the pinned do- 
main walls and also the possible role of the long-range 
dipolar interactions which exist in ferroelectric materi- 
als, it is thus necessary to perform a direct analysis 
of the static domain wall configuration, extracting the 
roughness exponent ^ and the effective domain wall di- 
mensionality dcs- Although measurements of this kind 
have been performed on other elastic disordered systems 
such as vortices (using neutron diffraction and decora- 
tion) [ll[l3, charge density waves (using X-ray diffrac- 
tion) lla]i contact lines [ly], and ferromagnetic domain 
walls [a] , a successful comparison between the experimen- 
tally observed roughness exponent and theoretical predic- 
tions could only be carried out in magnetic systems. In 
these systems, good agreement with the value ( = 2/3 
predicted for one-dimensional (line) domain walls in a 
random bond disorder was found 0. Quantitative stud- 
ies of these phenomena in other microscopic systems are 
therefore clearly needed. Epitaxial perovskite ferroelec- 
tric thin films with high crystalline quality and precisely 
controllable thickness are an excellent model system for 
such studies. 

In this paper we report on the first direct measurement 
of ferroelectric domain wall roughness. To image ferro- 
electric domains with the nanometer resolution required 
by such studies, atomic force microscopy (AFM) was used 
[ij, Il5[ l2fl I21I ] . Relaxation of the domain walls to their 
equilibrium configuration at short length scales allowed 
us to obtain values of 0.26 for the roughness exponent 
^. In the same films, the dynamic exponent /i was found 
to be ~ 0.5-0.6 from independent measurements of do- 
main wall creep. An analysis of these results gives an 
effective dimensionality of d ~ 2.5 for the domain walls, 
in good agreement with theoretical calculations for a two- 
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dimensional elastic interface in the presence of random- 
bond disorder and long range dipolar interactions (25j . 

The ferroelectric domain wall roughness studies were 
carried out in three c-axis oriented PbZro.2Tio.8O3 films, 
50, 66 and 91 nm thick, epitaxially grown on single crys- 
talline (001) NbiSrTiOs substrates by radio-frequency 
magnetron sputtering, as detailed in ref. H^H^- In these 
films, the polarization vector is parallel or anti-parallel to 
the c-axis and can be locally switched by the application 
of voltage signals via a metallic AFM tip, using the con- 
ductive substrate as a ground electrode l^^. The 
resulting ferroelectric domains are imaged by piezoforce 
microscopy (PFM) j2^. To measure domain wall rough- 
ness, we wrote linear domain structures with alternating 
polarization by applying alternating ±12 V signals while 
scanning the AFM tip in contact with the film surface. 
We chose line widths of 1 - 1.5 /im, and lengths of 8 - 
15 /im to ensure that domain wall images (2.5 x 2.5 and 
5x5 jim^) used in the study could be taken in the cen- 
tral regions, away from possible edge effects. Multiple 
domain structures were written in photolithographically 
pre-defined areas on each sample. More than 100 differ- 
ent ferroelectric domain walls were written and imaged 
in the three films. 

From these measurements, we extracted the correlation 
function of relative displacements 



B{L) = {[uiz + L) - u{z)]^) 



(1) 



where the displacement vector u{z) measures the defor- 
mation of the domain wall from an elastically optimal flat 
configuration due to pinning in favorable regions of the 
potential landscape. (• • •) and denote the thermody- 
namic and ensemble disorder averages, respectively. Ex- 
perimentally, the latter is realized by averaging over all 
pairs of points separated by the fixed distance L, rang- 
ing from 1 to 500 pixels (5 or 10 nm - 2.5 or 5.0 /im, 
depending on the image size) in our measurements. As 
shown in Fig. ^ for the three different films used, we ob- 
serve a power-law growth of B{L) at short length scales, 
comparable to the ~ 50-100 nm film thickness, followed 
by saturation of B{L) in the 100-1000 nm^ range [s^. 
The observed B{L) saturation indicates that the walls do 
not relax at large length scales from their initial straight 
configuration, dictated by the position of the AFM tip 
during writing. To ensure that domain wall relaxation 
was not hindered by the pinnin g pla nes of the lattice po- 
tential in the ferroelectric films 1^ , we wrote sets of 
domain walls at different orientations with respect to the 
crystalline axes in the 66 nm film. We found no corre- 
lation between the roughness of domain walls and their 
orientation in the crystal. This result is in agreement 
with previous studies 0| pointing out the negligible role 
of the commensurate potential compared to the effects of 
disorder. 

To investigate the possibility of thermal relaxation at 
ambient conditions, we then measured a set of domain 
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FIG. 1: Average displacement correlation function B{L) for 
different sets of ferroelectric domain walls in 50 (a), 66 (b) 
and 91 (c) nm thick films, shown out to L = 300 nm. Power 
law growth of B{L) is observed at short length scales, followed 
by saturation, suggesting a non-equilibrium configuration at 
large L. 



walls over a period of 1 month. No relaxation from the 
flat as-written configuration at large L is apparent visu- 
ally (Fig. Ufa)), or when comparing B{L) extracted for 
this set of domain walls at different times (Fig. I^b)). 
These data strongly indicate that ambient thermal acti- 
vation alone is not sufficient to equilibrate the domain 
walls over their entire length [s^. These results are in 
agreement with our previous studies [^Il5| in which both 
linear and nanoscopic circular domains remained com- 
pletely stable over 1-5 month observation periods. Such 
high stability is inherent to the physics of an elastic dis- 
ordered system, where energy barriers between different 
metastable states diverge as the electric field driving do- 
main wall motion goes to zero. This is an advantage 
for possible memory or novel filter applications 0, but 
also makes relaxation exceedingly slow. In fact, we be- 
lieve that even the relaxation leading to the observed 
power-law growth of B{L) at smaller length scales is not 
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FIG. 2: PFM images of the same set of domain walls taken 
(a) 1 hour after writing and (b) 1 week later. The white 
triangle indicates the same domain wall in each image, (c) 
The average B(L) for this set of domain walls. 



purely thermal, but occurs during the writing process it- 
self. When the direction of the applied electric field is 
reversed to form the alternating domain structure, the 
neighboring region already written with the opposite po- 
larity nonetheless experiences the resulting electric field, 
allowing the domain wall to locally reach an equilibrium 
configuration. 

From the power-law growth of B{L) at these short 
length scales, we extract a value for the roughness ex- 
ponent C- This exponent characterizes the roughness of 
the domain wall in the random manifold regime where an 
interface individually optimizes its energy with respect 
to the disorder potential landscape and B{L) scales as 
B{L) oc L^^. As shown in Fig. IHta), a linear fit of the 
lower part of the ln(_B(L)) vs ln(L) curve allows 2^ to be 
determined. Average values of C ~ 0.26, 0.29 and 0.22 
were obtained for the 50, 66 and 91 nm thick films, re- 
spectively, indicated by the horizontal lines in Fig.|3Ib). 

In addition to the investigations of static domain wall 
roughness described above, we independently measured 
domain wall dy namics in each film, using the approach 
detailed in 0, 0|. As shown in Fig. ^ we observe 
the non-linear velocity response to applied electric fields 
characteristic of a creep process, with values of 0.59, 0.58 
and 0.51 for the dynamical exponent /i in the 50, 66 and 
91 nm thick films, respectively [37l |. 

When these data are analyzed in the theoretical frame- 
work of a disordered elastic system, they provide infor- 
mation on the microscopic mechanism governing domain 
wall behavior. The direct measurement of domain wall 
roughness clearly rules out the lattice potential as a dom- 
inant source of pinning. In that case, the walls would 
have been fiat with B{L) ~ a^, where the lattice spac- 
ing a ~ 4 A is the period of the pinning potential 
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FIG. 3: (a) Typical In-ln plot of B{L). Fitting the linear part 
of the curve (left of the vertical line) gives 2(". (b) Average 
values of the characteristic roughness exponent ^ extracted 
from the equilibrium portion of the B{L) data are 0.26, 0.29 
and 0.22 in the 50, 66 and 91 nm thick samples, respectively, 
indicated by the horizontal lines in the figure. 
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FIG. 4: Domain wall speed as a function of the inverse applied 
electric field extracted from measurements of domain growth 
in the 50, 66 and 91 nm thick films. The data fit well to 
V ^ exp[— C/-B''] with = 0.59, 0.58 and 0.51, respectively. 



Given the stability and reproducibility of the wall po- 
sition over time, shown in Fig. |51 the effect of thermal 
relaxation on the observed increase of B{L) can also be 
ruled out. The measured roughness must thus be at- 
tributed to disorder. Two disorder universality classes 
exist, with different roughness exponents. Random bond 
disorder, corresponding to defects maintaining the sym- 
metry of the two polarization states would change only 
the local depth of the ferroelectric double well potential. 
Theoretically, this disorder would lead to a roughness ex- 
ponent Crb = 2/3 in deff = 1 and Crb ~ 0.2084(4 - d^s) 
for other dimensions. Random field disorder, correspond- 
ing to defects which locally asymmetrize the ferroelectric 
double well would favor one polarization state over the 
other. Such disorder would lead to a roughness exponent 
Crf = (4 — (ioff)/3 in all dimensions below four. Should 
the wall be described by standard (short range) elastic- 
ity, dcs in the above formulas is simply the dimension d 
of the domain wall (d = 1 for a line, d = 2 for a sheet). 
However, in ferroelectrics the stiffness of the domain walls 
and thus their elasticity under deformations in the direc- 
tion of polarization is different from the one for deforma- 
tions perpendicular to the direction of polarization be- 
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cause of long range dipolar interactions [25j . The elastic 
energy (expressed in reciprocal space) thus contains not 
only a short range term H — ^ J^q Cc\{q)u* iq)u{q) with 
cr^q^ but also a correction term due to the dipolar 
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interaction Cdp = + {^q^ + fg^) where y is 

the direction of the polarization. P is the ferroelectric 
polarization and e and eo are the relative and vacuum 
dielectric constants. Because qy now scales as qy ~ 
the effective dimensio n df, ff to use in the above formulas 
is des ~ d + 1/2 [2^, [sj. Using the above expressions 
for the roughness exponent we see that the measured 
C ~ 0.26 value would give deS > 3 for random field disor- 
der, ruling out this scenario. On the other hand random 
bond disorder would give doff 2.5 - 2.9, a much more 
satisfactory value, which is compatible with a scenario of 
two-dimensional walls (sheets) in random bond disorder 
with long range dipolar interactions. 

This conclusion can be independently verified by the 
dynamic measurements, since the creep exponent fi is 
related to the static roughness exponent ^ via ^ — 
The values of these two exponents from the 



independent static and dynamic measurements can there- 
fore be used to calculate dgg. For the 50, 66 and 91 nm 
thick films we find doff = 2.42, 2.49 and 2.47, respectively, 
in very good agreement with the expected theoretical 
value for a two-dimensional elastic interface in the pres- 
ence of disorder and dipolar interactions. Taken together, 
these two independent analyses provide strong evidence 
that the pinning in thin ferroelectric films is indeed due 
to disorder in the random bond universality class. The 
precise microscopic origin of such disorder is still to be de- 
termined. Preliminary studies of domain wall dynamics 
in pure PbTiOs films gave comparable results to those for 
PbZro.2Tio.8O3, suggesting that the presence of Zr in the 
solid solution is not a dominant factor, and that other de- 
fects presumably play a more significant role. Note that 
for the short-range domain wall relaxation observed, the 
walls are in the two-dimensional limit. However, if equi- 
librium domain wall roughness could be measured for 
larger L, a crossover to one-dimensional behavior would 
be expected, with a roughness exponent C = 2/3. 

In conclusion, we used AFM measurements of high 
quality epitaxial PbZro.2Tio.8O3 thin films to obtain the 
roughness exponent ( for ferroelectric domain walls. This 
is the first direct observation of static domain wall rough- 
ness in ferroelectric systems, and, combined with our 
measurements of domain wall dynamics, provides a co- 
herent physical image of their behavior in the framework 
of elastic disordered systems. 
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In other elastic disordered systems relaxation was 
achieved by the applying small oscillating driving forces 
[33, l3^ . or by subcritically driving the manifold into 
an equilibrium configuration |||. Using Pt electrodes de- 
posited on top of the writing areas we therefore applied 
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both DC and AC voltages to difFerent domain structures. 
Unfortunately, the large degree of leakage in such thin 
films, avoided in the local probe configuration, did not 
allow the domain walls to reach equilibrium at large L, 
with no relaxation observed when comparing B{L) for 



domain wall imaged before and after voltage application. 
[37] These ^ values are lower than the three values measured 
in 14], but consistent with all the subsequent measure- 
ments performed on nine other films, all grown under 
similar conditions. 



